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ABSTRACT 


Atmospheric  Wind  Measurements  Using 
a  50  MHz  Imaging  Doppler  Interferometer  Radar 


by 


Michael  Jerome  Volek,  Master  of  Science 
Utah  State  University 


Major  Professor:  Dr.  Gene  W.  Adams 
Department:  Soil  Science  and  Biometeorology 


I^The  50  megahertz  radar  located  near  LaSalle,  Colorado, 
also  known  as  MENTOR,  was  operated  as  an  Imaging  Doppler 
Interferometer  (IDI)  in  order  to  collect  radar  data.  These 
data  were  collected  from  18:25  Universal  Time  to  20:15 
Universal  Time  on  14  October  1988.  The  data  were  used  to 
measure  winds  in  the  lower  atmosphere  from  less  than  1 
kilometer  up  to  approximately  11.25  kilometers.  The  raw  data 
were  transformed  from  time  domain  to  frequency  domain  so  that 
the  IDI  technique  could  be  employed.  The  number  of  returns 
decreased  with  altitude,  as  did  their  radial  distance 


off-zenith.  Filters  were  developed  to  reject  points  that 
did  not  meet  certain  criteria,  and  then  algorithms  were 
developed  to  derive  wind  speeds  and  directions^  The  ^ 

horizontal  wind  and  its  zonal  and  meridional  components  were 
compared  to  rawinsonde  data  gathered  at  the  radar  site 


concurrently  with  the  IDI  measurements,  to  data  gathered  from 


ix 


the  ST  (stratosphere-troposphere)  radar  at  Platteville, 
Colorado,  and  to  the  14  Oct/12: 00  UT  and  15  0ct/00: 00  UT 
rawinsonde  runs  at  Denver.  The  vertical  velocities  measured 
were  compared  to  vertical  velocities  derived  from  the  Air 
Force  Global  Weather  Central  (AFGWC)  High  Resolution  (HIRAS) 
model  for  Denver.  Also,  baroclinicity  values  were  derived 
from  the  MENTOR  data  and  plotted  as  a  function  of  altitude. 

(101  pages) 


CHAPTER  I 


INTRODUCTION 

Scientists  have  long  been  interested  in  atmospheric 
winds  for  various  reasons.  At  first,  they  were  an  important 
consideration  in  weather  forecasting.  Today,  however,  the 
scope  of  importance  of  atmospheric  winds  is  much  broader  as 
man  utilizes  the  atmosphere  in  various  ways.  Commercial  and 
military  aviation,  the  Space  Shuttle,  ballistic  missiles,  and 
many  other  applications  depend  on  the  atmospheric  winds,  and 
the  ability  to  determine  the  atmospheric  wind  structure  is 
critical  to  the  success  of  these  systems.  Errors  in  wind 
measurement  can  have  tragic  and  costly  results.  For  example, 
commercial  airlines  need  wind  data  from  the  lower  atmosphere 
in  order  to  safely  and  economically  plan  their  routes.  That 
is  why  many  new  methods  of  instrumentation  are  being 
developed — to  help  meet  the  growing  need  for  atmospheric  wind 
data. 

Man  has  used  many  different  methods  to  sample  the  upper 
air  winds.  The  first  was  by  using  a  rawinsonde,  which 
involves  attaching  a  telemetering  system  to  a  balloon,  then 
using  simple  trigonometry  to  calculate  the  winds  (Byers, 
1974)  .  The  upper  limit  of  the  rawinsonde  is  3  6  km.  Some 
innovative  design  changes  in  the  balloon  resulted  in  the 
Jimsphere,  a  balloon  that  is  overfilled  so  that  it  expands 
on  ascent.  Jimspheres  can  remain  aloft  at  fixed  altitudes 
and  telemeter  wind  information  to  ground  sites.  Their  upper 
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limit  is  27  km.  After  their  development,  rockets  were  used 
to  gather  wind  data  from  regions  inaccessible  to  balloons. 
Scientists  could  now  sample  the  upper  atmosphere  up  to  240 
km  (Byers,  1974)  . 

Radar  was  a  technological  breakthrough  that  was  of  great 
benefit  to  atmospheric  scientists.  Incoherent  scatter  radars 
measure  the  plasma  bulk  flow  of  the  upper  atmosphere  by 
sensing  weakly  re-radiated  energy  from  upper  atmospheric  free 
electrons.  The  radar  then  measures  the  frequency  change  due 
to  the  Doppler  shift,  and  then  this  frequency  change  can  be 
used  to  measure  radial  velocities.  Various  techniques  can 
then  be  employed  to  determine  the  horizontal  and  vertical 
velocities.  This  type  of  radar  can  only  receive  returns  from 
ionospheric  heights  (above  60  km)  .  The  mesosphere- 
stratosphere-troposphere  radar  (MST)  uses  a  large,  coherent, 
monostatic  radar  operating  in  the  VHF-UHF  range  to  obtain 
echoes  in  the  optically  clear  atmosphere  from  1-100  km 
(Balsley,  1981) .  The  primary  echoing  mechanism  involves 
scattering  from  turbulent  irregularities  having  a  scale  size 
equal  to  one-half  the  radar  wavelength.  MST  radars  transmit 
three  beams — one  vertically  and  two  off-zenith.  The  radar 
then  receives  echoes  from  refractive  index  variations  in  the 
scattering  volume  being  scanned.  The  drawback  to  this 
technique  is  that  it  depends  on  upon  the  existence  of  these 
variations  at  the  three  transmit  angles  simultaneously. 
Since  this  is  not  always  the  case,  MST  radars  must  sometimes 
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sample  a  region  over  a  long  period  of  time.  The  MST  radar 
measures  the  Doppler  shift  induced  by  the  moving  scatterers 
to  determine  winds  (Balsley,  1981) . 

Spaced  Antenna  Drift  (SAD)  radars  are  incoherent 
systems.  They  also  measure  Doppler  shifts  induced  by 
refractive  index  variations.  The  transmitting  radar  beam  is 
directed  vertically,  and  three  or  more  spaced  receiving 
antennas,  also  directed  vertically,  are  used  to  observe  the 
return  echoes.  The  fluctuations  of  echo  strength  at  the 
spaced  antennas  are  recorded.  In  the  presence  of  a 
horizontal  wind,  these  fluctuations  show  relative  time 
displacements.  Horizontal  winds  can  be  computed  from  the 
time  displacements  (Briggs,  1980) .  This  method  does  not 
depend  on  the  simultaneous  existence  of  refractive  index 
variations,  and  is  the  technique  on  which  the  IDI  is  based. 

Meteor  radars  detect  the  ionized  trail  left  by  a 
decaying  meteor  as  it  disintegrates  in  the  upper  atmosphere 
above  80  km.  This  limits  how  low  in  the  atmosphere  these 
radars  can  measure  winds.  Meteor  echoes  are  manifested  as 
events  with  very  sharp  increases  in  amplitude  followed  by 
exponential  decays.  The  phase  history  is  well  defined  for 
the  echo  lifetime,  and  easily  yields  the  radial  Doppler 
velocity  of  the  echo.  One  can  calculate  the  mean  wind 
profile  in  a  given  altitude  interval  if  there  are  three  or 
more  echoes  present  in  the  interval  (Brosnahan  et  al.,  1987) . 

The  fact  that  MST  radars  may  require  long  sampling  times 
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makes  them  susceptible  to  the  influence  of  long  period  wave 
structures  (such  as  gravity  waves) ,  and  mesoscale  motions. 
This  makes  the  wind  profiles  obtained  by  them  questionable 
(Royrvik,  1983) .  However,  comparisons  of  tropospheric  and 
stratospheric  winds  determined  by  MST  radars  and  rawinsondes 
show  good  agreement  (e.g.  Warnock  et  al.,  1978  and  Fukao  et 
al.,  1982)  .  MST  radar  measurements  of  the  mesosphere  did  not 
compare  as  well  to  rocketsonde  data,  even  though  there  was 
fair  correlation  between  the  methods  (Smith  and  Fritts, 
1984) .  SAD  radar  measurements  of  the  troposphere  and 
stratosphere  also  show  good  agreement  with  rawinsonde  data 
(Rottger  and  Ierkic,  1985) . 

Each  of  these  methods  has  provided  results  in  obtaining 
wind  measurements  for  certain  atmospheric  regions.  Figure 
1  shows  the  various  techniques  used  in  measuring  winds  in  the 
atmosphere  and  the  altitudes  from  which  each  technique 
provides  results  (Gage  and  VanZandt,  1981) . 

A  new  technique  has  been  developed  using  a  50  MHz  radar 
as  an  Imaging  Doppler  Interferometer.  The  frequency  is  high 
enough  to  allow  clear-air  returns  from  the  stratosphere  and 
troposphere  and  can  be  used  to  measure  winds  in  these 
atmospheric  regions.  A  big  advantage  to  this  system  is  that 
it  does  not  require  the  scatterers  to  be  in  a  certain  place 
at  a  certain  time.  This  radar  illuminates  a  large  volume  of 
scatterers,  and  can  translate  the  return  signal  into 
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parameters  one  can  use  to  calculate  winds.  Long  sampling 
times  are  not  required  to  obtain  a  large  number  of  echoes, 
thus  making  the  measurements  more  reliable.  This  system 
consists  of  thirty-two  Yagi  antennas;  one  array  of  sixteen 
is  used  as  a  transmitter  and  the  other  array  of  sixteen  is 
used  as  a  receiver.  This  method  uses  Fourier  transforms  to 
convert  time  domain  data  into  the  frequency  domain.  One  can 
then  use  the  frequency  information  to  obtain  the  parameters 
which  are  used  to  calculate  wind  profiles  (Adams  et  al., 
1986)  . 

The  purposes  of  this  thesis  are  as  follows: 

Investigate  the  use  of  a  50  MHz  Imaging  Doppler 
Interferometer  (IDI)  radar  to  measure  tropospheric 
winds. 

-  Determine  the  best  combination  for  weighting  and/or 
rejecting  scattering  points. 

Plot  the  resulting  wind  profiles  with  height  and  compare 
them  with  data  from  the  Denver  rawinsonde  runs,  the  ST 
radar  located  near  Platteville,  CO,  and  the  rawinsonde 
runs  made  concurrently  at  the  radar  site.  I  will  also 
investigate  any  discrepancies  between  the  results 
obtained  by  IDI  and  the  other  methods. 

If  the  results  are  inconclusive,  I  will  report  such  a 
finding,  and  suggest  further  research  into  the  IDI 
technique . 
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CHAPTER  II 

LITERATURE  REVIEW 

Probing  of  the  atmosphere  by  radar  depends  on  the 
scattering  and  reflection  from  spatial  refractive  index 
variations  at  scales  of  one-half  the  radar  wavelength  (Gage 
and  Balsley,  1980) .  Depending  on  the  altitude  under 
consideration,  these  variations  can  correspond  to  changes  in 
pressure,  temperature,  humidity,  and  electron  density.  These 
changes  are  what  causes  a  radar  signal  to  be  scattered  or 
reflected  (Rottger,  1984) .  In  the  troposphere,  electron 
density  can  be  neglected,  and  in  the  upper  troposphere, 
humidity  is  small  enough  to  be  negligible  (Gage  et  al., 
1981) .  These  changes  in  the  refractive  index  can  best  be 
illustrated  by  the  following  equation  (Gage  and  Balsley, 
1980)  : 

n  -  1  =  (3.73  x  10_1e)/T2  +  (7.76  x  10'6P)/T  -  Ne/2NC  (1) 
where  n  is  the  refractive  index,  P  is  the  atmospheric 
pressure  in  millibars,  e  is  the  partial  pressure  of  water 
vapor  in  millibars,  T  is  the  absolute  temperature  in  ’K,  N„ 
is  the  electron  number  density  in  m'3,  Nc  is  the  critical 
plasma  density  in  m'3  for  the  radar  frequency,  where  Nc  =  1.24 
x  10  2f 2  and  f  is  the  radar  frequency  in  Hz. 

Doppler  radar  works  in  the  following  manner:  when  a 
pulse  of  electromagnetic  radiation  impinges  on  a  target,  it 
forces  molecular  vibrations  in  synchronism  with  the 
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time-changing  electric  and  magnetic  fields.  If  the  target 
is  stationary,  its  molecules  will  vibrate  at  the  frequency 
of  the  radiation  field.  If  the  target  is  moving  toward  the 
transmitter  at  a  velocity  v,  its  vibrational  frequency  is 
higher  by  v/A  ,  (where  A  is  wavelength)  because  the  target 
molecules  experience  more  rapid  fluctuations  of  electric  and 
magnetic  force  due  to  the  target  experiencing  more  frequent 
strikes  by  the  radar  pulse  when  it's  moving  toward  the  radar 
(Battan,  1973)  .  This  is  the  Doppler  effect.  The  vibrating 
molecules  in  turn  generate  electromagnetic  fields,  which  then 
radiate  outward  from  the  target.  For  a  monostatic  radar 
where  the  transmitter  and  receiver  are  collocated,  the 
frequency  of  the  scattered  radiation  is  Doppler-shifted;  its 
Doppler  frequency  is  given  by  fd  =  -2v*/A,  where  vR  is  the 
radial  (line-of-sight)  component  of  the  velocity,  positive 
being  away  from  the  radar.  Similar  reasoning  applies  for  a 
target  moving  away  from  the  transmitter  due  to  the  Doppler 
effect  (Doviak  and  Zrnic,  1984) . 

Radar  techniques  have  proven  to  be  a  very  powerful  tool 
for  studying  the  earth's  atmosphere.  A  technique  that  uses 
the  scattering  property  of  water/ ice  particles  and/or  the 
fluctuations  of  the  index  of  refraction  in  clear  air  have 
been  used  to  study  the  mesosphere  and  stratosphere  (Woodman 
and  Guillen,  1974) .  The  echoes  obtained  from  these  regions 
come  from  backscattering  produced  by  fluctuations  of  the 
atmospheric  dielectric  constant.  Analysis  of  the  returned 
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signals  yields  information  on  the  dynamics  of  the  large-scale 
motion  as  well  information  on  the  degree  of  turbulence 
(Woodman  and  Guillen,  1974) . 

The  dielectric  properties  and,  therefore,  the 
corresponding  fluctuations,  are  produced  by  different 
phenomena  within  the  troposphere,  stratosphere  and  the 
mesosphere.  In  the  troposphere  and  stratosphere,  the 
dielectric  properties  are  determined  by  the  density  (i.e., 
by  its  temperature  at  a  given  pressure) ,  while  in  the 
mesosphere  these  properties  are  determined  by  the  number  of 
free  electrons.  However,  the  results  obtained  by  radar 
probing  are  independent  of  the  mechanisms  responsible  for  the 
fluctuations.  If  a  dielectric  medium  is  illuminated  by  an 
electromagnetic  wave,  the  wave  propagates  according  to  the 
refractive  and  attenuating  properties  of  the  medium.  In  any 
dielectric,  no  matter  how  transparent,  a  fraction  of  the 
power  carried  in  the  primary  wave  is  scattered  in  all 
directions.  Fluctuations  in  a  medium  are  a  random  process 
in  space  and  time  (Woodman  and  Guillen,  1974) . 

Early  experiments  postulated  that  clear  air  turbulence 
in  the  troposphere  was  responsible  for  the  enhancement  of 
dielectric  fluctuations.  Turbulence  in  a  region  with  a 
temperature  inversion  is  very  efficient  in  producing 
fluctuations  in  temperature,  and  hence  in  the  index  of 
refraction.  The  vertical  gradient  of  the  index  of  refraction 
is  given  by  (Woodman  and  Guillen,  1974) : 
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dN/dZ  =  -79  X  10'6p/T2  X  (dT/dz  +  7.)  (2) 

where  p  is  the  pressure  in  millibars,  T  is  the  absolute 
temperature,  and  7,  is  the  adiabatic  lapse  rate  in  "K/m. 
Under  most  circumstances,  the  scattering  is  due  to 
atmospheric  refractive  index  variations  arising  from  locally 
homogeneous,  isotropic  turbulence.  In  the  lower  troposphere, 
these  variations  can  also  be  caused  by  regular  variations  of 
both  temperature  and  humidity,  and  are  isotropic  in  nature 
(Rottger  and  Liu,  1978) .  Fresnel  refection  is  also 
responsible  for  producing  atmospheric  returns.  These  are 
caused  by  sharp  coherent  gradients  in  the  atmospheric 
refractive  index  that  extend  horizontally  and  are  produced 
by  stratified  clear  air  refractivity  structures.  These 

structures  are  usually  found  in  the  upper  troposphere  and  the 
stratosphere  (Balsley  and  Gage,  1981) .  Anisotropic 
turbulence  may  also  play  a  part  in  the  reflection  of  radar 
signals  (Rottger  et  al.,  1981). 

Pertinent  scale  sizes  of  the  turbulence  structures  which 
radar  can  detect  are  those  corresponding  to  one-half  the 
transmitted  wavelength.  Turbulent  flow  energy  is  transferred 
according  to  the  Kolmogorov  hypothesis,  which  says  that 
energy  enters  the  turbulence  spectrum  at  large  wavelengths 
and  is  transferred  to  smaller  wavelengths  until  it's 
dissipated.  This  is  called  energy  cascade  (Vinnichenko  et 
al.,  1980).  So,  the  largest  scale  sizes  correspond  to  the 
region  where  turbulence  is  generated  due  to  mean  flow 


11 


instability.  This  is  the  energy-containing  subrange. 
Wavelengths  are  on  the  order  of  several  kilometers  (Panofsky 
and  Dutton,  1984) .  The  buoyancy  subrange  follows,  which  is 
a  sink  for  turbulent  energy.  It  has  a  characteristic  scale 
size  of  several  hundred  meters  to  one  kilometer  (Vinnichenko 
et  al.,  1980).  Since  these  scales  are  larger  than  the 
wavelength  of  VHF  radars,  they  do  not  contribute  to  the 
scattering/reflection  process  (Rottger,  1984) .  Energy  from 
the  energy-containing  and  buoyancy  subranges  is  transmitted 
through  the  inertial  subrange.  Here,  there  is  no  gain  or 
loss  of  energy.  Wavelengths  are  on  the  order  of  several  to 
tens  of  meters  (Panofsky  and  Dutton,  1984).  It's  in  the 
inertial  subrange  where  VHF  radar  echoes  are  received 
(Rottger,  1984)  .  The  short  wavelength  limit  of  the  inertial 
subrange  is  determined  by  viscous  damping  and  is  given  by  the 
Kolmogorov  microscale  Ak: 

Ak  =  2ir  (i/3/ e) 1/4  (3) 

where  e  is  the  eddy  dissipation  rate  of  turbulence  and  v  is 
kinematic  viscosity  (Balsley,  1978) .  Scale  sizes  less  than 
the  microscale  comprise  the  viscous  subrange,  where  energy 
is  transferred  from  the  inertial  subrange.  This  energy  is 
converted  into  heat  due  to  viscosity  (Vinnichenko  et  al., 
1980) .  Wavelengths  in  this  subrange  are  on  the  order  of  1 
cm  in  the  troposphere,  and  on  the  order  of  1  m  in  the 
mesosphere.  Wavelengths  smaller  than  the  microscale  are 
strongly  damped,  so  that  the  microscale  essentially 
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determines  the  upper  frequency  limit  for  a  radar  system 
designed  to  investigate  a  specific  atmospheric  region 
(Balsley,  1978) .  However,  the  microscale  increases  with 
height,  since  it  is  inversely  proportional  to  the  density. 
Therefore,  the  wavelength  of  the  turbulence  structures  in 
the  inertial  subrange  increases  with  height.  This  is  why  MF 
radars  can  detect  echoes  in  the  mesosphere  (Rottger,  1984)  . 

The  return  signal  strength  in  a  region  decreases  with 
increasing  height  up  to  about  50  km,  then  begins  to  increase. 
Below  50  km,  the  decrease  with  height  is  primarily  due  to  the 
decreasing  atmospheric  density  and  the  increased  stability; 
above  50  km,  the  scattering  is  due  partly  to  the  increase 
with  height  of  the  ambient  ionization  near  the  bottom  of  the 
ionosphere  and  partly  due  to  enhanced  neutral  turbulence 
within  the  mesosphere  due  to  increased  viscous  effects 
(Balsley,  1978)  .  To  sum  it  all  up,  since  scattering 
typically  derives  from  turbulence  with  scale  sizes  equal  to 
one-half  the  transmitted  wavelength,  it's  necessary  to  have 
turbulent  structure  of  this  dimension  present  in  the 
scattering  volume  in  order  to  obtain  an  echo  (Balsley,  1981)  . 

Clear  air  radars  normally  require  MF-UHF  (300  kHz-3  GHz) 
transmission  frequencies,  peak  power  of  a  few  megawatts,  and 
antenna  apertures  of  hundreds  or  thousands  of  square  meters. 
The  systems  must  be  frequency  coherent  (capable  of 
determining  the  phase  of  the  incoming  signal  relative  to  the 
transmitted  pulse)  in  order  to  measure  the  Doppler 


13 


characteristics  of  the  echo.  One  must  assume  that  the 
scattering  volume  is  completely  filled  with  scatterers,  and 
that  the  medium  is  homogeneous  and  isotropic.  This 
assumption  holds  well  in  the  troposphere  and  mesosphere, 
where  homogeneous,  isotropic  turbulence  plays  a  major  role 
in  scattering  radar  pulses.  In  the  stratosphere,  however, 
the  stability  and  horizontal  stratification  of  refractive 
index  gradients  cause  the  returns  to  be  specular.  In  either 
case,  there  must  be  a  way  for  the  system  to  pick  out  real 
returns  from  noise.  There  are  various  ways  in  which  to 
detect  a  signal  from  the  atmosphere  in  the  presence  of  noise. 
Figure  2  shows  an  example  of  this.  The  returned  signal  is 
denoted  by  the  enhanced  peak  at  Sr(Af)  +  S„(Af)  ,  where  Sr  and 
SN  are  the  power  spectral  densities  of  the  signal  and  noise, 
respectively.  The  noise  power  spectral  density  in  the 
absence  of  a  signal  is  denoted  by  SN (Af )  ,  and  the  mean  value 
of  the  noise  is  indicated  by  the  horizontal  line  SN.  The 
noise  fluctuation  level  aSn  is  also  shown.  The  ±AF  limits 
are  determined  by  the  transmitted  pulse  rate  and  the  number 
of  coherent  integrations,  which  involves  the  summing  of  a 
number  of  samples  obtained  from  the  echoing  volume  after 
sequential  transmitter  pulses  and  prior  to  performing  a 
spectral  analysis  (Balsley,  1978) . 
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Fig.  2.  Plot  of  spectral  power  density  vs.  Doppler 
shift.  (From  Balsley,  1978.) 
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To  determine  the  vector  wind  at  a  given  height,  at  least 
three  Doppler  measurements,  each  in  a  different  direction, 
are  required.  Normally,  one  beam  is  directed  to  the  north, 
5 ° -15 0  off-zenith,  one  is  directed  to  the  east,  also  5“-15° 
off-zenith,  and  one  is  directed  vertically  (Rottger,  1984) . 
Sometimes  the  vertical  beam  is  omitted  on  the  assumption  that 
vertical  velocities  may  be  negligible;  this  will  provide  a 
2-dimensional  horizontal  velocity  (Balsley,  1981) .  The  two 
or  three  regions  from  which  Doppler  returns  are  received  are 
separated  from  each  other  in  space  and  sometimes  also  in  time 
(Waterman  et  al.,  1985).  Recently,  it  was  discovered  that 
there  seems  to  exist  short  period  fluctuations  with  spatial 
scales  less  than  the  distance  between  the  beam  positions. 
These  fluctuations  can  deteriorate  wind  velocity  estimations. 
Further  error  would  also  be  caused  by  the  wind  nonuniformity 
within  the  range  volume  which  is  determined  by  antenna 
aperture,  transmitted  pulse  width,  and  range  gate  sampling 
(Fukao  et  al.,  1988).  The  measured  velocity  is  assumed  to 
be  an  average  of  wind  velocity  within  the  range  volume,  where 
radio  wave  scattering  occurs.  The  average  does  not  reflect 
the  mean  wind  velocity  at  a  specific  height,  but  is  weighted 
by  turbulence  intensity  (Fukao  et  al.,  1988).  Wind 
estimation  error  is  then  caused  by  the  finiteness  of  the 
range  volume,  called  the  finite  range  volume  effect  (Fukao 
et  al.,  1988).  This  effect  can  produce  false  vertical  wind 
shears  of  the  horizontal  velocity. 
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Fresnel  reflection  from  a  stably  stratified  atmosphere 
is  considered  as  one  mechanism  which  causes  signal  power 
decreases  with  increased  zenith  angle.  This  is  known  as 
aspect  sensitivity  (Rottger,  1984) .  Another  mechanism 
involves  diffuse  reflection  from  a  rough  boundary  between 
laminar  and  turbulent  regions.  One  final  mechanism  involves 
the  incoherent  addition  of  reflected  signals  from  many  layers 
filling  the  volume  illuminated  by  the  radar  beam  (Fresnel 
scattering)  (Gage  et  al.,  1981;  Tsuda  et  al.,  1986). 

It  has  been  shown  that  the  aspect  sensitivity  of  a 
reflecting  region  increases  with  increasing  stability  (Tsuda 
et  al.,  1986).  This  can  account  for  the  high  aspect 
sensitivity  (specularity)  noted  in  the  stratosphere.  This 
means  that  the  stratosphere  is  filled  with  aspect  sensitive 
echoes.  Isotropic  echoes  independent  of  zenith  angle  are 
detected  only  in  the  troposphere.  Rottger  et  al.  (1981) 
determined  that  the  angular  dependence  for  aspect  sensitive 
upper  tropospheric  and  stratospheric  scatterers  was  about 
1.5-2. 5  dB/degree.  The  stratospheric  echoes  with  large 
aspect  sensitivity  cannot  be  attributed  to  the  contribution 
from  a  single  isolated  reflecting  surface  caused  by  an  abrupt 
vertical  change  in  the  refractive  index,  but  to  contributions 
from  many  reflecting  layers  that  fill  the  whole  stratosphere 
(Tsuda  et  al.,  1986).  Therefore,  measurements  at  large 
zenith  angles  should  be  avoided  in  order  to  minimize  errors 
in  the  horizontal  wind  velocities  (Tsuda  et  al.,  1986). 
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The  radar  interferometer  exploits  the  Doppler 
information  inherent  in  coherent  spaced  antenna  measurements. 
Cross  spectral  analysis  of  the  signals  received  at  separate 
antennas  can  be  used  to  locate  and  track  irregularities  while 
Doppler  sorting  can  be  used  to  discriminate  against  multiple 
targets,  provided  they  have  different  radial  velocities 
(Vincent,  1984) . 

The  interferometer  technique  allows  one  to  measure 
simultaneously  the  angular  spectrum  of  the  returns.  It  gives 
useful  information  on  wave  and  turbulence  structure,  and 
improves  the  vertical  velocity  measurements  (Rottger  and 
Ierkic,  1985).  The  technique  works  only  for  small,  localized 
scatterers  or  reflectors,  and  not  for  horizontally  stratified 
clear  air  refractivity  structures. 

To  get  a  value  of  the  vertical  velocity,  one  must  know 
the  inclination  of  the  reflecting  structures.  The 
interferometer  technique  allows  one  to  obtain  the  most 
accurate  inclination  measurements  as  well  as  the  most 
accurate  aspect  sensitivity  measurements  (Rottger  and  Ierkic, 
1985) .  The  average  inclination  of  the  reflecting  structures 
can  give  an  estimate  of  the  inclination  of  isentropic 
surfaces  (Rottger,  1984) .  Large  scale  vertical  motion  is  a 
consequence  of  nearly  horizontal  motion  on  isentropic 
surfaces,  since  parcels  of  air  in  adiabatic  motion  conserve 
potential  temperature  (Gage,  1983) .  If  one  can  measure  the 
motion  along  an  isentropic  surface,  then  the  velocity  vector 
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can  be  broken  down  into  vertical  and  horizontal  components. 
If  this  3-dimensional  wind  vector  is  denoted  by  VM/  then, 

| VM|  =  (VH2  +  w2  )1/2,  where  VH  is  the  horizontal  velocity 
component  and  w  is  the  vertical  velocity  component.  The 
horizontal  velocity  can  be  broken  down  into  zonal  (u)  and 
meridional  (v)  components,  or,  VH  =  (u2  +  v2)1/2.  Once  the 
horizontal  and  vertical  components  are  known,  one  can  compute 
the  baroclinicity  S  by  the  simple  trigonometric  relation  <5 
=  tan’1  (w/VH)  .  Baroclinicity  is  used  by  meteorologists  to 
determine  baroclinic  instability,  which  effects  synoptic 
scale  disturbances  that  travel  along  the  jet  stream  and 
develop  into  weather-producing  systems  (Holton,  1972) .  Radar 
measurements  of  baroclinicity  could  be  used  by  numerical 
weather  prediction  models  as  an  initial  condition  of  the 
atmosphere.  Radar  interferometry  provides  the  best 
measurement  of  vertical  velocity,  and  could  therefore  be  used 
to  determine  baroclinicity  (Ierkic  and  Rottger,  1984) .  The 
vertical  velocity  spectra  measured  are  due  to  waves,  but  the 
horizontal  velocity  spectra  have  been  attributed  to  both 
quasi-two-dimensional  turbulence  and  to  the  spectrum  of 
buoyancy  waves  (Gage  and  Nastrom,  1984) . 

Imaging  Doppler  Interferometry 

The  technique  used  to  analyze  the  data  presented  herein 
involve  the  principles  of  imaging  Doppler  interferometry. 
A  50  MHz  radar  operating  as  an  Imaging  Doppler  Interferometer 
(IDI)  was  used,  and  the  general  concepts  of  interferometry 
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are  borrowed  from  optics.  The  IDI  depends  on  atmospheric 
motions  to  Doppler-shift  the  individual  scatterers  so  that 
they  can  be  separated  later  by  performing  a  Fourier 
transform. 

Radar  interferometry  involves  making  use  of  the  phase 
difference  between  signals  received  at  two  or  more  spatial 
locations.  At  radio  frequencies  used  to  measure  atmospheric 
properties,  one  can  measure  the  phases  of  the  received 
voltages  on  two  or  more  antennas  simultaneously.  High-speed 
processing  then  allows  one  to  make  conventional  radar 
measurements . 

The  IDI  separates  multiple  targets  by  first  transforming 
the  data  from  time  domain  to  frequency  domain,  then  using  the 
differences  in  Doppler  phase  spectra  to  locate  the  individual 
targets.  The  time  series  of  complex  voltages  from  each 
antenna  are  Fourier  transformed,  and  then  individual  spectral 
features  are  spatially  located  from  the  phase  difference 
between  antennas.  The  complex  Fourier  transform  returns  both 
an  amplitude  and  phase  at  each  Doppler  frequency. 

Figure  3  illustrates  how  a  single  target,  moving  at  a 
horizontal  velocity  VH,  is  measured.  The  location  and 
velocity  of  the  target  can  be  measured  by  Fourier- 
transforming  the  time  series  of  the  two  complex  voltages  to 
get  two  power  and  phase  spectra.  Figure  4  is  an  example  of 
the  power  and  phase  spectra  for  the  moving  target  in  Figure 
3.  Both  power  spectra  identify  the  same  Doppler  frequency 
f,  which  can  be  used  to  obtain  the  the  radial  velocity  VR: 
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Fig.  4.  Complex  Fourier  spectra  for  the  two  antennas 
in  Figure  3,  showing  the  signature  of  a  single  moving 
target.  (From  Adams,  et  al.,  1985.) 
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VR  =  f  c/2F  =  VHsin0  (4) 

where  f  is  the  Doppler  frequency,  F  is  the  radar  frequency, 
and  c  is  the  speed  of  light.  The  difference  in  phase  angles 
can  be  used  to  calculate  the  zenith  angle: 

e  =  sin'1  (0A/2jt1)  (5) 

where  1  is  the  spacing  between  the  two  antennas,  A  is  the 
radar  wavelength,  and  <p  is  the  phase  (Adams  et  al.,  1986). 

The  50  MHz  IDI,  however,  uses  thirty-two  antennas,  so 
a  slightly  different  approach  is  used.  Radar  returns  are 
received  on  eight  independent  antennas.  Each  antenna's  data 
are  separately  Fourier-transformed,  and  the  resulting  complex 
Fourier  spectra  at  each  range  gate  are  examined  for  Doppler 
bins  whose  sources  can  be  located  by  spatial  interferometry. 
Real  sources  are  identified  by  requiring  agreement  among  the 
independent  antenna  pairs  on  the  two-dimensional  angular 
location  of  each  source. 

One  can  then  determine  a  set  of  characteristics  that 
will  define  each  individual  scattering  point;  these 
characteristics  are  the  scattering  point  parameters  (SPPs) 
and  consist  of: 

z  :  the  altitude  at  which  the  echo  is  located; 

VR  :  the  echo  radial  velocity; 

dm  :  the  echo-location  angle  in  the  E-W  plane; 

0NS  :  the  echo-location  angle  in  the  N-S  plane; 

Vgy,  :  the  average  voltage  on  the  E-W  antenna; 

VNS  :  the  average  voltage  on  the  N-S  antenna; 
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<Pm  :  the  phase  on  the  E-W  antenna,  relative  to  the 
center  of  the  array; 

0NS:  the  phase  on  the  N-S  antenna,  relative  to  the 
center  of  the  array. 

Now  we  have  defined  a  three-dimensional  coordinate 
system,  which  one  can  transform  to  a  Cartesian  system  (X  is 
east,  Y  is  north,  and  Z  is  vertical) .  The  scattering  point 
parameters  are  sorted  by  altitude.  One  can  now  calculate  the 
mean  apparent  motion  of  each  scattering  point.  The  radial 
velocity  of  each  scatterer  can  be  calculated  from  equation 
4  as: 

Vrj  =  cf  j/2F  (6) 

where  j  =  1,2,3, . ,J,  where  J  is  the  number  of  scattering 

points.  One  can  also  write  the  vector  radial  velocity  as: 

Vrj  =  VrjIrj  (7) 

where  1RJ  is  a  unit  vector  in  the  radial  direction  passing 
through  the  jth  scattering  point,  given  by: 

lpj  =  ljlx  +  fijlj  (8) 

where  lj,  mj,  and  nj  are  the  direction  cosines  of  the  jth 
scattering  point,  given  by: 

1  =  sin  9m  (9) 

m  =  sin  (10) 

n  =  (l-l2  -mJ  )uz  (11) 

Then,  the  mean  apparent  motion  vector  can  be  denoted  as: 

Vm  =  ulx  +  vly  +  wl2  (12) 

Three  scattering  points  are  the  minimum  required  to  determine 
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the  three  components  of  motion.  With  these  points,  one  can 
calculate  a  least-squares  fit  to  the  SPPs  as  follows:  Since 

vRj  =  '  Irj  (13) 

for  all  j  if  the  rms  fit  is  perfect,  the  rms  error  is  given 
by: 

e  =  {1/JECVrj  -  ( ul j  +  vrrtj  +  wnj)  ] 2  }1/2  (14) 

We  need  to  chose  the  values  of  u,  v,  and  w  that  minimize  e, 
or: 


de/du  =  de/dv  -  de/dw  =  0  (15) 

which  then  yields: 

uElj2  +  vSl^  +  wSl^nj  =  ZVrjIj  (16) 

uSljmj  +  vSitij2  +  wZmjnj  =  EV^m.,  (17) 

uZljmj  +  vZmjnj  +  wZn^  =  ZVRjnj  (18) 


These  equations  can  then  be  solved  simultaneously  for 
u,  v,  and  w,  which  then  give  the  three-dimensional  bulk  flow 
for  each  range  gate. 
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CHAPTER  III 

DATA  ANALYSIS 

The  MENTOR  50  MHz  IDI  Radar 

The  data  used  for  the  analysis  were  collected  by  the  50 
MHz  MENTOR  Imaging  Doppler  Interferometer  (IDI)  radar, 
located  near  LaSalle,  Colorado  (40 ° 16'  42"N,  104 ° 35'  55"W) . 

The  MENTOR  (Meteor  Echoes?  No  Transmitter,  Only  Receivers) 
had  originally  been  designed  to  be  a  passive  system, 
receiving  pulses  reflected  out  of  ST  radar  beams  by  meteor 
trails.  A  50  kw  transmitter  was  later  added  so  that  a 
stand-alone  IDI  radar,  used  to  measure  tropospheric  and 
stratospheric  winds,  could  be  deployed.  The  MENTOR  operates 
at  frequencies  between  49.8  and  49.92  MHz,  with  50  kw 
peak-pulse-power  in  pulse  lengths  of  1.67  or  6.67  /isec.  The 
pulse  repetition  rate  is  variable  from  470  juHz  to  12  kHz. 
One  complete  sounding  contains  256  samples,  and  it  takes  10 
seconds  to  obtain  one  sounding. 

Figure  5  shows  graphically  how  the  system  is  set  up. 
Two  identical  arrays  of  16  5-element  vertically  directed 
Yagis  are  used;  one  array  is  the  transmitter,  the  other  is 
the  receiver.  Both  are  arrayed  in  a  four-by-four  matrix  with 
1.05-0.7-1.05  wavelength  spacing  between  both  the  row  and 
column  antennas.  The  transmit  antennas  are  fed  with  a  16-way 
power  splitter  with  equal  power  and  phase.  The  resulting 
beam  is  approximately  20°  wide  at  the  3-dB  points.  On  the 
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receive  array,  the  16  signals  are  power-divided  with  0° 
hybrids,  and  these  signals  are  recombined  with  4 -way  power 
combiners  to  obtain  a  total  of  8  outputs  consisting  of  4  row 
arrays  and  4  column  arrays.  Each  antenna's  pattern  is 
fan-shaped,  being  20“  wide  along  the  line  of  Yagis  and  about 
50°  perpendicular  to  the  line  of  Yagis.  The  MENTOR  uses  a 
direct-conversion  receiver  technique.  A  total  of  8  channels 
are  used  to  amplify  and  detect  the  8  antenna  signals 
(Brosnahan  et  al.,  1988). 

During  data  collection,  coherent  averaging  is  done  and 
is  dynamically  adjusted  to  maintain  a  TDA  output  rate  of  256 
samples  every  10  seconds.  An  FFT  is  then  done  on  the  data 
to  yield  results  in  the  frequency  domain,  from  which  the 
scattering  point  parameters  are  derived.  From  this  point, 
one  can  then  calculate  wind  profiles  in  the  lower  and  middle 
atmosphere. 

Tropospheric  Wind  Analysis 

In  the  past,  a  2.66  MHz  IDI  radar  was  used  to  determine 
mesospheric  wind  profiles.  The  MENTOR  was  designed  to 
measure  winds  from  the  troposphere  and  stratosphere  as  well. 
It  can  measure  winds  in  these  regions  due  to  its  high 
sensitivity  and  its  ability  to  gather  the  data  needed  to 
provide  wind  calculations  in  a  relatively  short  time  period. 
The  MENTOR  also  possesses  extremely  fast  data-processing 
hardware,  which  makes  gathering  the  data  and  its  conversion 
to  scattering  point  parameters  via  an  FFT  routine  possible. 
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Figure  6  is  a  computer  block  diagram  of  the  MENTOR  system 
setup. 

In  order  to  calculate  the  winds,  I  first  developed 
software  to  make  the  basic  calculations.  I  tested  it  in  two 
different  ways.  First,  I  took  three  points  at  5.1  km.  One 
point  was  directly  overhead  ( 6^  and  0NS  both  equal  to  zero), 
one  was  10°  off-zenith  on  the  E-W  array  and  the  other  was 
10°  off-zenith  on  the  N-S  array.  The  test  wind  was  30  m/s 
from  260°.  This  gave  a  zonal  component  (u)  of  29.54  m/s  and 
a  meridional  component  (v)  of  5.21  m/s.  I  let  the  vertical 
component  (w)  equal  zero  for  this  test.  I  then  derived 
radial  velocities  for  each  of  the  three  points,  and  used 
these  values  along  with  the  values  of  and  0NS  in  the 
program.  I  obtained  z  =  5.1  km,  u  =  29.6  m/s,  v  =  5.24  m/s, 
w  =  0.0  m/s,  horizontal  velocity  =  30.06  m/s,  and  horizontal 
direction  =  260°.  This  provided  a  simple  check  to  insure  the 
logic  and  correctness  of  the  routines.  I  then  analyzed 
synthetic  data  derived  to  further  test  this  program  (see 
Adams  et  al.  (1985)  for  information  relating  to  the 
generation  of  these  data) .  The  horizontal  velocity  was  input 
to  be  a  constant  30  m/s,  with  the  horizontal  direction  vector 
rotating  counterclockwise. 

In  order  to  analyze  the  synthetic  data,  I  first 
calculated  the  rms  error  using  the  technique  described  in  the 
last  chapter  (equation  14)  .  After  calculating  the  components 


Fig.  6.  MENTOR  block  system  diagram. 
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and  the  rms  error  for  each  range  gate,  I  then  discarded  those 
points  where  the  individual  error  was  greater  than  twice  the 
rms  error  for  the  particular  range  gate  and  recalculated  the 
winds.  The  input  to  the  synthetic  data  (Figure  7)  is  the 
desired  result.  I  varied  the  number  of  points  per  velocity 
in  each  profile  to  see  if  more  points  produced  a  better 
average.  The  results  for  10,  15,  30,  and  60  points  are  shown 
in  Figures  8-11.  Each  profile  approximates  the  input  quite 
well.  It  appears  that  more  points  per  velocity  calculation 
provides  better  average  results,  with  the  curve  becoming 
smoother  with  each  subsequent  increase  in  points.  The 
largest  value  of  VH  for  10  points  was  32.5  m/s,  for  60  points 
it  was  31.8  m/s;  the  smallest  value  of  VH  for  10  points  was 
27.2  m/s,  for  60  points,  it  was  29.1  m/s.  However,  it  is 
apparent  that  these  differences  are  minor,  and  that  using 
fewer  points  provides  better  altitude  resolution.  The  ideal 
situation  would  be  to  use  as  few  points  as  possible  so  that 
one  can  obtain  a  very  dense  profile. 

The  actual  data  analyzed  were  obtained  14  Oct  88  from 
18:25  UT  to  20:15  UT.  These  data  will  be  compared  to 
National  Weather  Service  rawinsonde  runs  at  14  0ct/12:00  UT 
and  15  Oct/00: 00  UT  out  of  Denver,  two  rawinsonde  balloons 
launched  at  the  radar  site  to  coincide  with  the  radar  data 
acquisition,  and  data  from  the  ST  radar  run  by  NOAA/WPL  near 
Platteville,  CO.  The  rawinsonde  launches  at  the  radar  site 
are  extremely  desirable:  Denver  is  approximately  60  miles  SW 
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Fig.  7.  Wind  component  input  to  synthetic  data 
generator. 
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Fig.  8.  Test  calculations  using  synthetic  data  and  10 
points  per  velocity. 
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Fig.  9.  Test  calculations  using  synthetic  data  and  15 
points  per  velocity. 
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Fig.  10.  Test  calculations  using  synthetic  data  and  30 
points  per  velocity. 
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Fig.  ll.  Test  calculations  using  synthetic  data  and  60 
points  per  velocity. 
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of  the  radar  site,  in  the  foothills  of  the  Rockies,  and  the 
winds  measured  there  can  be  influenced  more  by  the  mountains 
than  those  at  the  radar  site,  which  lies  on  the  high  plateau 
east  of  the  base  of  the  mountains.  Also,  winds  are  a  dynamic 
feature;  even  5  minutes  can  make  a  substantial  difference  in 
wind  velocity.  Thus  any  discrepancies  noted  between  the 
Denver  runs  and  the  IDI  could  be  cleared  up  by  the  data  from 
the  local  rawinsonde  run.  A  comparison  between  the  IDI  and 
the  conventional  ST  radar  will  also  prove  interesting. 
Platteville  is  approximately  10  miles  SW  of  LaSalle,  and  it's 
also  on  the  high  plains  east  of  the  foothills.  There  should 
be  little  difference  in  the  mean  flow  over  the  two  sites; 
therefore  this  can  be  considered  a  comparison  of  the  two 
radar  techniques.  In  comparing  the  IDI  data  with  the  LaSalle 
rawinsonde  data,  one  must  consider  the  fact  that  as  the 
rawinsonde  balloon  rises,  it  is  carried  horizontally  with  the 
wind.  On  the  average,  a  balloon  takes  one  hour  to  reach  a 
height  of  5  km.  If  one  uses  an  average  horizontal  wind  of 
5  m/sec  over  that  altitude  range,  it  is  displaced  by  43  km 
from  its  launch  site.  Then,  the  wind  measured  at  5  km  by  the 
IDI  is  overhead,  while  that  measured  by  the  balloon  may  now 
be  43  km  away. 

The  two  balloon  launches  at  LaSalle  were  at  18:22  and 
19:42  UT,  and  the  Platteville  ST  run  was  at  19:00  UT.  I 
developed  software  to  break  down  the  horizontal  velocities 
from  the  ST  and  balloon  observations  into  components;  these 
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components  will  be  compared  to  those  obtained  by  MENTOR. 

Dr.  Gene  Adams  (Personal  Communication,  1988)  developed 
the  main  software  program  that  I  used  to  derive  the 
tropospheric  winds  from  MENTOR.  One  of  the  major  problems 
one  needs  to  resolve  is  determining  which  points  received  are 
good  and  eliminating  any  ambiguities  in  the  reception  of 
points.  For  this  purpose,  the  antennas  were  set  as  shown  in 
Figure  5.  The  two  Yagis  closest  to  each  other  are  0.7 A 
apart,  while  the  distance  between  the  end  and  middle  Yagis 
are  1.05A,  making  the  total  distance  between  the  end  Yagis 
2 . 8 A .  Thus,  one  can  measure  the  phase  at  each  of  the  four 
Yagis:  <pt  and  04  at  the  ends  and  <p2  and  03  at  the  two  in  the 
middle.  There  are  also  four  complex  voltages  measured,  which 
gives  three  independent  phase  differences:  0x-02  (1.05A), 

03— 02  ( 0 . 7  A ) ,  and  04-03  (1.05A).  Each  phase  difference 

corresponds  to  one  or  more  possible  zenith  angles  through  the 
interferometer  equation: 

A  0  =  -27rDsin0  (19) 

where  D  is  the  antenna  spacing  in  wavelengths.  Thus,  for  the 
two  pair  of  Yagis  1.05A  apart,  A0=2 . 17rsin0 ,  and  for  the  pair 
that  are  0.7A  apart,  A0=1. 47rsin0 .  Solutions  of  these  two 
equations  yield  several  values  of  5.  A  variable  zenith  angle 
window  is  then  used.  There  are  three  independent 
measurements  of  the  E-W  zenith  angle,  which  will  spread 
across  some  range  of  zenith  angles,  and  three  simultaneous 
measurements  of  the  N-S  zenith  angle  and  their  resultant 
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spread.  The  sum  of  the  two  maximum  spreads  of  zenith  angles 
in  each  direction  derived  from  the  three  independent  phase 
differences  (the  0.7A  pair  and  the  two  1.05A  pair)  must  be 
less  than  the  specified  zenith  angle  window  for  one  of  the 
possible  ±2 n  permutations.  The  basic  design  of  the  window 
is  as  follows:  there  are  several  zenith  angles  in  both  the 
E-W  and  N-S  directions.  For  the  first  time  through  the 
window  routine  (the  E-W  values) ,  the  maximum  allowed  zenith 
angle  spread  is  the  user-defined  zenith  angle  window.  The 
second  time  through  the  routine  (the  N-S  values) ,  the  maximum 
allowed  spread  is  the  user-defined  window  minus  the  actual 
E-W  value.  Then,  equation  19  gives  the  relationship  between 
the  antenna-to-antenna  phase  difference  and  zenith  angle. 
For  spacings  larger  than  0.5A,  two  or  more  zenith  angles  are 
possible  solutions  for  a  given  phase  difference.  For  the 
1.05A  and  the  0.7A  spacings,  three  solutions  exist.  One 
determines  the  maximum  and  minimum  zenith  angles  for  each 
direction,  and  calculates  the  spread,  or: 

^spread  —  ®  max  ~  ^min  (20) 

If  the  smallest  spread  is  greater  than  the  zenith  angle 
window,  the  point  is  rejected.  If'  it's  less,  then  the  actual 
zenith  angle  has  been  located.  The  spreads  calculated  on 
both  the  E-W  and  N-S  arrays  are  written  as  scattering  point 
parameters.  A  consensus  angle  is  then  defined  as  the  middle 
of  the  calculated  window,  or: 

e 


consensus 


(*m.*  -  *ain)/2 


(21) 
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A  second  check  is  then  done  by  coherently  averaging  the 
voltages  from  each  pair  of  antennas  located  1.05 A  apart,  then 
calculating  the  phase  difference  between  them.  This  will 
generate  five  values  of  $  from  a  single  phase  difference. 
The  values  of  6  obtained  by  both  methods  are  then  compared; 
only  one  of  the  solutions  obtained  by  coherent  averaging  will 
be  within  the  window  centered  on  the  consensus  zenith  angle. 
This  is  the  final  zenith  angle,  which  improves  the  accuracy 
of  the  calculated  consensus  value  by  a  factor  of  (NJ  1/2*NC*D, 
where  is  the  number  of  incoherent  averages,  Nc  is  the 
number  of  coherent  averages,  and  D  is  the  antenna  spacing  in 
wavelengths.  If  there  are  no  values  that  lie  within  the 
window  centered  on  the  consensus  angle,  then  the  point  is 
rejected.  Setting  the  window  to  a  larger  value  obviously 
allows  more  noise  into  the  system,  but  also  lets  more  true 
points  in  as  well.  Smaller  values  decrease  both  noise  and 
true  points.  One  must  use  trial-and-error  methods  to 
determine  which  value  of  the  window  is  optimum,  and  also  make 
use  of  filters  to  weed  out  system-generated  points  and  any 
others  that  may  contaminate  the  data.  For  this  analysis,  I 
found  that  a  zenith  angle  window  of  24*  worked  best. 

I  used  a  linear  polarization  filter  in  conjunction  with 
doing  an  FFT  on  the  time  domain  data.  This  is  done  to  match 
the  transmit  polarization. 

After  the  time-domain  data  were  transformed  to  the  SPP 
format,  I  merged  consecutive  strings  of  10-second  soundings 
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together  to  create  one  long  sounding  of  various  lengths.  I 
did  this  so  that  I  could  plot  the  various  parameters  against 
each  other  to  see  if  there  were  any  patterns  evident  in  the 
distribution.  I  also  wanted  to  have  various  sounding  lengths 
available  when  I  did  the  wind  profiling;  this  way  I  could 
determine  what  the  optimum  sounding  length  should  be.  The 
sounding  lengths  were  2,  3,  4,  6,  8,  and  10  minutes.  Initial 
review  of  10  minutes  of  data  taken  from  18:25:10  to  18:35:10 
UT  on  14  Oct  88  showed  that  the  radial  velocities  obtained 
by  MENTOR  were  contaminated  with  noise.  Figure  12  shows  a 
scatter  plot  of  radial  velocity  versus  radial  zenith  angle, 
showing  radial  velocities  much  greater  than  zero  for  angles 
near  zenith.  This  would  translate  roughly  into  high  vertical 
velocities,  on  the  order  of  ±30-40  m/s,  since  the  radial 
velocity  near  zenith  is  mostly  composed  of  vertical  velocity. 
Such  high  vertical  velocities  are  unrealistic,  and  would 
correspond  to  horizontal  velocities  near  zero  for  points  near 
zenith.  In  order  to  filter  out  these  points,  six  10-second 
soundings  were  coherently  averaged,  and  the  results  plotted. 
A  least-squares  polynomial  was  then  fit  to  these  data.  This 
filtered  out  the  high  radial  velocities  close  to  zenith  while 
keeping  the  other  points.  The  coherent  average  was  done  in 
order  to  enhance  the  pattern,  especially  near  zenith.  The 
one-minute  average  cut  down  the  number  of  points  and  made  it 
easier  to  determine  a  pattern  and  to  fit  the  curve.  The 
algorithm  was  then  used  as  a  filter  for  the  merged  data 
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Scatter  plot  of  radial  zenith  angle  vs.  radial 
from  10  minutes  of  MENTOR  data. 
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files. 

Generally,  points  that  are  received  near  zenith  give  a 
good  representation  of  vertical  velocities,  while  points 
further  off-zenith  will  provide  good  results  for  horizontal 
velocities.  Knowing  this,  I  separated  the  merged  data  into 
bins;  the  first  one  for  points  that  lie  within  ±2°  of  zenith, 
another  for  points  between  ±2°  and  ±10°  off-zenith,  and 
another  for  points  between  ±10°  and  ±25°  off-zenith.  I 
decided  to  try  using  only  points  within  ±2°  of  zenith  for 
vertical  velocity  profiles,  and  points  further  off-zenith  for 
horizontal  profiles.  Interestingly,  I  obtained  fairly  good 
results  for  both  vertical  and  horizontal  velocities  using  the 
data  within  ±2°  of  zenith,  while  the  results  using  points 
further  off-zenith  provided  results  that  did  not  correlate 
at  all  with  the  results  from  the  other  techniques.  Once  I 
noted  this  phenomenon,  I  tried  opening  up  the  point- 
collection  window  to  ±3°,  since  there  would  be  more  points 
in  that  extra  degree  that  I  could  use  to  generate  more  wind 
values  in  a  profile.  However,  this  only  contaminated  the 
data,  and  I  found  no  agreement.  I  then  opened  up  the  window 
to  ±10°  using  a  1°  increment.  The  results  were  the  same:  no 
agreement.  I  then  tried  keeping  only  points  within  ±1°  of 
zenith,  but  there  were  not  enough  points  to  provide  more  than 
2  wind  values  in  the  troposphere.  At  first,  it  seemed  very 
unusual  that  only  those  points  within  ±2°  of  zenith  would 
provide  reasonable  wind  profiles. 
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I  plotted  the  east-west  and  north-south  zenith  angles 
versus  altitude  and  the  east-west  and  north-  south  returned 
voltages  versus  the  east-west  and  north-south  zenith  angles, 
respectively  (Figures  13-16) .  These  scatter  plots  show  what 
is  happening:  note  the  cluster  of  scattering  points  near 
zenith  in  Figures  13  and  14.  This  shows  that  a  majority  of 
returns  are  coming  from  near  zenith.  The  scatter  plots  in 
Figures  15  and  16  show  that  the  received  voltages  are  orders 
of  magnitude  higher  close  to  zenith.  This  means  that  all  the 
returns  from  points  not  close  to  zenith  are  noise.  This  is 
due  to  the  fact  that  the  power  transmitted  falls  off  rapidly 
with  zenith  angle,  so  that  there  is  effectively  no  signal 
going  out  past  approximately  ±2  ° ,  and  many  of  the  returns 
past  that  are  noise.  Figures  17  and  18,  which  are  plots  of 
the  east-west  and  north-south  received  voltages  versus 
altitude,  further  illustrate  this. 

From  Figures  13  and  14,  one  notices  that  there  are  not 
as  many  scattering  points  near  zenith  at  altitudes  between 
4  and  10  km  as  there  are  at  altitudes  below  4  km  and  above 
10  km.  Then  Figures  17  and  18  show  that  the  returned 
voltages  have  peaks  below  4  km  and  above  10  km.  This  can  be 
attributed  to  the  loss  of  signal  power  at  points  far  enough 
off-zenith.  This  signal  power  fall-off  is  a  function  of  the 
antenna  frequency;  the  higher  the  frequency,  the  narrower  the 
transmitted  beam.  One  could  lower  the  frequency  and  widen 
the  beam,  but  at  the  expense  of  low-altitude  returns.  The 
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Fig.  14.  Scatter  plot  of  N-S  zenith  angle  vs.  altitude 
from  MENTOR. 
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Fig.  15.  Scatter  plot  of  E-W  received  voltage  vs.  E-W 
zenith  angle  from  MENTOR. 
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Fig.  18.  Scatter  plot  of  N-S  received  voltage  vs 
altitude  from  MENTOR. 
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occurrance  of  good  returns  only  near  zenith  could  be  a 
property  of  the  scatterers,  independent  of  beamwidth. 
Incidentally,  other  experiments  using  radar  interferometers 
to  measure  atmospheric  winds  found  that  even  at  small 
off-zenith  angles,  the  horizontal  component  of  the  average 
velocity  is  not  negligible  and  could  be  measured  (Rottger  and 
Ierkic,  1985) . 

After  calculating  multiple  wind  profiles  and  doing 
comparisons  with  the  other  data,  I  was  able  to  determine  that 
using  a  zenith  angle  window  of  12°  or  16°  did  not  provide 
enough  data  to  plot  a  usable  wind  profile,  even  after  merging 
enough  consecutive  files  to  create  a  10-minute  sounding. 
Opening  the  window  to  20°  and  24°  lets  in  more  points  (and 
noise) ,  but  the  filters  that  are  applied  to  the  data  work  to 
effectively  keep  the  noise  out.  When  calculating  the 
profiles  using  these  data,  I  discovered  that  using  fewer 
points  allows  for  better  altitude  resolution,  so  more  points 
can  be  plotted  with  height.  I  also  noted  that  the  results 
I  obtained  by  discarding  points  that  had  an  error  of  more 
than  twice  the  rms  error  in  the  given  range  gate  were  not 
significantly  different  than  those  obtained  by  keeping  all 
the  points.  Even  after  opening  the  zenith  angle  window  to 
24°,  the  filtering  process  takes  its  toll  on  the  number  of 
points  left  with  which  to  calculate  wind  profiles.  For  this 
reason,  I  used  10  minutes  of  consecutive  data  to  calculate 
the  profiles  presented  in  the  next  chapter.  There  are  simply 
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not  enough  points  when  using  soundings  much  shorter  than 
that.  I  used  two  separate  10-minute  soundings  to  calculate 
profiles,  the  first  from  18:25:10  to  18:35:10  UT  and  the 
second  from  18:45:10  to  18:55:10  UT.  These  were  the  only  two 
tapes  in  which  there  were  enough  continuous  data  to  complete 
a  10-minute  sounding. 

If  more  continuous  data  were  available,  longer  soundings 
could  be  created,  and  profiles  with  better  resolution  could 
be  calculated.  However,  this  too  has  its  drawbacks.  As  a 
scatterer  approaches  the  radar  beam  from  a  distance,  it  moves 
through  several  Fourier  windows  as  it  moves  across  the  sky 
and  through  the  beam.  This  means  that  the  Doppler  shift 
changes  each  time  the  scatterer  passes  through  a  different 
window.  The  longer  the  sounding,  the  more  windows  the  target 
passes  through.  This  has  an  effect  on  the  wind  velocity  that 
is  calculated,  and  is  known  as  velocity  smearing. 
Unfortunately,  smearing  is  a  complex  process,  and  no  one  to 
date  has  been  able  to  determine  how  to  "desmear"  the  data. 

After  the  data  analysis  was  completed,  I  found  out  that 
there  was  a  problem  with  one  of  the  antennas.  Apparently, 
the  y-component  on  antenna  number  four  was  always  zero,  which 
would  induce  errors  in  the  echo  location.  Even  though  this 
is  a  problem,  so  much  filtering  takes  place  that  most  of  the 
bad  pints  were  probably  removed.  One  could  also  develop  a 
filter  to  eliminate  all  the  data  from  antenna  four  before 
using  them  to  calculate  winds. 
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CHAPTER  IV 

RESULTS  AND  DISCUSSION 

This  chapter  contains  the  wind  profiles  calculated  using 
the  MENTOR  data  and  the  technique  developed  to  analyze  them. 
The  wind  velocities  calculated  are  not  instantaneous,  but 
time-averaged  bulk  flow  rates.  Only  extremely  sensitive 
sonic  anemometers  can  measure  instantaneous  wind  velocities. 
The  MENTOR  is  providing  the  average  wind  velocity  in  a  given 
range  gate  over  a  10  minute  period. 

The  data  presented  include  plots  of  the  horizontal 
velocity,  its  zonal  and  meridional  components,  vertical 
velocity,  the  direction  from  which  the  horizontal  wind  is 
coming,  and  baroclinicity .  Standard  meteorological  notation 
is  used  in  order  to  coincide  with  the  data  taken  from  the 
rawinsondes  and  the  ST  radar.  This  means  that  the  horizontal 
direction  of  the  wind,  a,  is  given  by  a  =  tan'1  (u/v) ,  where 
u  is  the  zonal  component  and  v  is  the  meridional  component. 
A  wind  blowing  straight  from  the  north  is  0°,  from  the  east 
is  90”,  from  the  south  is  180”,  and  from  the  west  is  270”. 
The  vertical  velocity  is  either  positive  or  negative,  with 
positive  being  upwards  and  negative  being  downwards.  This 
convention  is  adopted  assuming  the  bulk  flow  is  along 
isentropic  surfaces.  These  surfaces  generally  slope  upwards 
in  the  Northern  Hemisphere,  so  a  bulk  flow  up  a  surface  is 
ascending,  and  is  associated  with  instability,  while  flow 
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down  a  surface  is  descending,  and  is  associated  with  stable 
conditions  (Gage,  1983) . 

The  baroclinicity  <5  is  also  positive  or  negative. 
Positive  <5  is  associated  with  instability,  while  negative  <5 
is  associated  with  stability.  It  is  a  relative  measure  of 
the  ratio  of  vertical  to  horizontal  velocity,  with  a  larger 
positive  S  corresponding  to  greater  instability  and  a  larger 
negative  S  corresponding  to  greater  stability. 

The  plots  of  the  horizontal  speed  and  direction,  and  the 
zonal  and  meridional  components  can  be  easily  compared  to  the 
LaSalle  and  Denver  rawinsonde  runs  and  the  Platteville  ST 
radar  run.  However,  this  is  not  the  case  with  the  vertical 
velocity,  since  none  of  these  methods  measures  vertical 
velocity.  This  is  because  vertical  velocity  is  difficult  to 
measure.  ST  radars  do  not  usually  employ  a  vertical  beam  to 
measure  this  component,  since  the  magnitudes  are  generally 
quite  small,  and  can  be  assumed  to  be  zero  in  many  cases 
v  ..age,  1983)  .  Since  there  are  no  direct  measurements  of 
vertical  velocities  for  comparison  with  MENTOR  values,  I  used 
data  derived  by  the  U.  S.  Air  Force  Global  Weather  Central's 
High  Resolution  (HIRAS)  Model  for  Denver  at  18:00  UT.  This 
data  is  derived  from  other  fields  measured  by  the  12:00  UT 
Denver  rawinsonde,  which  is  then  input  into  a  numerical  model 
to  generate  the  values  for  18:00  UT.  There  are  also  no 
baroclinicity  measurements  to  use  as  a  comparison;  as  a 
matter  of  fact,  I  have  not  found  a  plot  of  baroclinicity 
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anywhere  in  the  literature. 

To  generate  these  profiles,  I  used  soundings  comprised 
of  10  minutes  of  consecutive  data.  The  zenith  angle  window 
was  set  at  24°  to  allow  in  as  many  points  as  possible.  I 
used  the  linear  polarization  and  least-squares  polynomial 
filters  mentioned  in  the  last  chapter,  and  used  points  within 
±2°  of  zenith.  The  profiles  generated  don't  contain  as  many 
points  as  the  rawinsonde  or  ST  radar  runs,  but  do  contain 
enough  to  paint  a  general  picture  of  the  atmospheric 
structure  below  12  km.  Also,  the  horizontal  winds  measured 
by  the  rawinsondes  are  normally  accurate  to  within  ±1  m/s 
(Warnock  et  al.,  1978),  while  there  is  no  mention  in  the 
literature  of  the  accuracy  of  the  ST  data.  The  accuracy  of 
the  rawinsonde  method  was  challenged  by  Fukao  and  his 
co-workers  (Fukao  et  al.,  1982),  who  argued  that  differences 
between  winds  measured  by  MST  radar  and  rawinsondes  at 
Arecibo,  Puerto  Rico,  were  due  mainly  to  inherent 
experimental  errors  in  the  rawinsonde  technique.  They  went 
on  to  mathematically  prove  their  theory. 

Figures  19-22  show  comparisons  of  the  horizontal  speed 
and  direction,  and  zonal  and  meridional  components  from 
MENTOR  data  taken  between  18:25:10  and  18:35:10  UT  with  the 
LaSalle  rawinsondes  and  the  Platteville  ST  radar.  Note  the 
terms  LORAN  and  OMEGA;  these  refer  to  two  different 
techniques  used  to  determine  wind  velocities  using 
rawinsondes.  There  is  good  agreement  between  the  profiles 
of  horizontal  speed  and  direction.  One  must  remember  that 
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Fig.  19.  A  comparison  of  horizontal  wind  speeds 
determined  by  the  18:25:10-18:35:10  MENTOR  data,  the 
LaSalle  rawinsondes,  and  the  Platteville  ST  radar. 
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Fig.  20.  A  comparison  of  horizontal  directions 
determined  by  the  18:25:10-18:35:10  MENTOR  data,  the 
LaSalle  rawinsondes,  and  the  Platteville  ST  radar. 
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Fig.  21.  A  comparison  of  zonal  wind  speeds  determined 
by  the  18:25:10-18:35:10  MENTOR  data,  the  LaSalle 
rawinsondes,  and  the  Platteville  ST  radar. 
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Fig.  22.  A  comparison  of  meridional  wind  speeds 
determined  by  the  18:25:10-18:35:10  MENTOR  data,  the 
LaSalle  rawinsondes,  and  the  Platteville  ST  radar. 
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the  resolution  on  the  MENTOR  data  is  not  as  good  as  it  could 
be  if  the  system  was  able  to  measure  points  further 
off-zenith.  The  MENTOR  profile  follows  the  general  trend 
given  by  the  other  techniques.  There  is  a  large  discrepancy 
in  direction  above  10  km,  but  one  must  always  keep  in  mind 
that  each  of  these  soundings  was  taken  at  a  different  time 
using  different  techniques  and  types  of  equipment.  Not  only 
are  there  experimental  errors  which  are  inherent  in  any 
measurement,  but  the  atmosphere  is  extremely  dynamic. 
Turbulent  eddies  are  created  and  destroyed  in  time  scales  on 
the  order  of  seconds.  These  eddies  transfer  horizontal 
momentum  vertically,  which  produces  fluctuations  in  the 
horizontal  wind  flow.  Therefore  one  can  expect  to  see 
differences  in  the  profiles;  the  rawinsonde  and  ST  data  do 
not  even  agree  totally.  The  discrepancy  in  the  direction  is 
due  to  the  difference  in  the  zonal  component  measured  by 
MENTOR.  It  measured  a  small  easterly  component  (speeds  less 
than  zero  are  easterly)  where  the  others  measured  relatively 
large  westerlies.  It's  interesting  to  note  the  good 
agreement  in  measurements  of  the  meridional  component  in 
Figure  22.  The  strong  directional  shear  noted  in  Figure  20 
and  the  corresponding  shear  in  the  meridional  component  could 
be  an  indication  of  an  upward-propagating  mountain  lee  wave. 
Figures  23-26  show  the  same  MENTOR  data  compared  against  the 
Denver  rawinsonde  data.  Once  again,  there  is  fairly  good 
agreement  between  the  profiles,  even  though  6  hours  separate 
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Fig.  23.  A  comparison  of  horizontal  wind  speeds 
determined  by  the  18:25:10-18:35:10  MENTOR  data  and  the 
Denver  rawinsondes. 
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Fig.  24.  A  comparison  of  horizontal  directions 
determined  by  the  18:25:10-18:35:10  MENTOR  data  and  the 
Denver  rawinsondes. 
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Fig.  25.  A  comparison  of  zonal  wind  speeds  determined 
by  the  18:25:10-18:35:10  MENTOR  data  and  the  Denver 
rawinsondes . 
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Fig.  26.  A  comparison  of  meridional  wind  speeds 
determined  by  the  18:25:10-18:35:10  MENTOR  data  and  the 
Denver  rawinsondes. 
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the  MENTOR  profile  from  each  Denver  one.  Also,  the  two  sites 
are  about  60  miles  apart.  Denver  is  at  the  base  of  the 
Rockies  while  LaSalle  is  on  the  high  plains.  There  are  more 
discrepancies  in  all  the  profiles,  but  that  should  be 
expected  in  this  case.  Figure  27  is  a  comparison  of  vertical 
velocity  from  MENTOR  with  those  derived  from  the  HIRAS  model. 
These  are  very  good  results,  since  they  both  show  upward 
vertical  motion  and  the  magnitudes  are  similar.  This  result 
would  indicate  that  the  MENTOR  is  well  suited  to  measure 
vertical  velocities.  Figure  28  shows  a  plot  of  baroclinicity 
from  MENTOR.  The  high  positive  value  near  2  km  would 
indicate  strong  instability,  with  increasing  stability  above 
that  until  just  below  10  km. 

Figures  29-32  are  comparisons  of  the  wind  profiles 
generated  by  MENTOR  data  from  18:45:10  to  18:55:10  UT,  the 
LaSalle  rawinsondes  and  the  Platteville  ST  radar.  Since 
there  were  not  as  many  points  in  this  10  minute  sounding  as 
there  were  in  the  previous  one,  there  are  not  as  many  points 
in  its  profile.  This  may  also  explain  why  the  results  are 
not  as  good  as  those  obtained  10  minutes  earlier.  For 
example,  Figure  29  shows  a  large  increase  in  winds  between 
3.6  km  and  5.2  km.  It's  possible  that  this  would  not  appear 
if  there  were  more  data  to  work  with;  more  data  would  improve 
the  resolution.  However,  the  data  does  follow  the  general 
trend  well,  especially  the  increase  at  5.2  km,  which  was  also 
detected  by  the  LORAN  rawinsonde  out  of  LaSalle.  The 
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Fig.  27.  A  comparison  of  vertical  velocities  determined 
by  the  MENTOR  18:25:10-18:35:10  data  and  the  AFGWC 
HIRAS  model  for  Denver. 
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Fig.  29.  A  comparison  of  horizontal  wind  speeds 
determined  by  the  MENTOR  18:45:10-18:55:10  data,  the 
LaSalle  rawinsondes,  and  the  Platteville  ST  radar. 
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Fig.  30.  A  comparison  of  horizontal  directions 
determined  by  the  MENTOR  18:45:10-18:55:10  data,  the 
LaSalle  rawinsondes,  and  the  Platteville  ST  radar. 
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Fig.  31.  A  comparison  of  zonal  wind  speeds  determined 
by  the  MENTOR  18:45:10-18:55:10  data,  the  LaSalle 
rawinsondes,  and  the  Platteville  ST  radar. 
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Fig.  32.  A  comparison  of  meridional  wind  speeds 
determined  by  the  MENTOR  18:45:10-18:55:10  data,  the 
LaSalle  rawinsondes,  and  the  Platteville  ST  radar. 
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horizontal  directions  obtained  this  time  around  to  not  agree 
as  well  as  those  10  minutes  earlier.  This  is  mainly  caused 
by  the  discrepancy  in  meridional  components  measured  by 
MENTOR  and  by  the  other  techniques  (see  Figure  32)  .  The  only 
good  agreement  is  at  the  lowest  and  highest  altitudes. 
Otherwise,  MENTOR  measures  high  southerly  components,  while 
the  others  measure  relatively  small  northerly  components. 
MENTOR  and  the  LORAN  rawinsonde  do  correlate  somewhat  on  the 
strong  southerly  component  at  5.2  km.  Figure  30  shows  that 
the  directions  measured  by  MENTOR  are  for  the  most  part 
within  90°  of  those  measured  by  the  other  methods.  Figure 
31  shows  that  there  is  good  agreement  on  the  zonal  components 
measured  by  each  method.  The  same  arguments  used  in 
describing  the  first  set  of  MENTOR  data  comparisons  can  be 
used  here  as  well.  It  demonstrates  how  turbulent  and  chaotic 
the  atmosphere  is.  Figures  33-36  show  this  MENTOR  profile 
compared  to  those  from  the  Denver  rawinsondes.  Figure  33 
shows  some  discrepancy  in  the  horizontal  wind  measured  at  5.2 
km,  but  otherwise  good  agreement.  Once  again,  better 
resolution  could  solve  this  problem.  The  horizontal 
directions  in  Figure  34  once  again  show  the  MENTOR  profile 
offset  from  the  others  by  about  30 '-60°;  note  the  strong 
directional  shear  measured  by  the  12:00  UT  rawinsonde.  This 
is  an  indication  of  how  quickly  the  atmosphere  can  change, 
and  is  possibly  another  indication  of  an  upward-propagating 
mountain  lee  wave.  It  is  possible  that  the  MENTOR  does  pick 
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Fig.  33.  A  comparison  of  horizontal  wind  speeds 
determined  by  the  MENTOR  18:45:10-18:55:10  data  and  the 
Denver  rawinsondes. 
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Fig.  34.  A  comparison  of  horizontal  directions 
determined  by  the  MENTOR  18:45:10-18:55:10  data  and  the 
Denver  rawinsondes. 
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Fig.  35.  A  comparison  of  zonal  wind  speeds  determined 
by  the  MENTOR  18:45:10-18:55:10  data  and  the  Denver 
rawinsondes . 
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Fig.  36.  A  comparison  of  meridional  wind  speeds 
determined  by  the  MENTOR  18:45:10-18:55:10  data  and  the 
Denver  rawinsondes. 
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up  this  wave  motion,  but  due  to  the  long  sounding  times  and 
the  effects  of  smearing,  does  not  show  the  sharp  distinction, 
but  rather  a  smoother  directionalvariation  with  height.  One 
can  see  from  Figures  3  5  and  3  6  that  there  is  pretty  good 
agreement  between  the  zonal  components  but  not  as  much 
between  the  meridional  components.  Figure  37  shows  the 
comparison  between  vertical  velocities  from  MENTOR  and  the 
HIRAS  model.  There  is  once  again  good  agreement  here,  and 
it  appears  that  the  MENTOR  has  detected  a  wave  with  a 
vertical  wavelength  of  about  8  km.  One  could  make  a  case  for 
this  being  a  gravity  wave  that  was  generated  by  the  Rockies, 
since  not  only  is  there  vertical  evidence  of  propagation  but 
possibly  some  evidence  of  horizontal  propagation  as  well. 
Figure  38  is  a  plot  of  baroclinicity ,  and  shows  that  the 
atmosphere  is  stable  between  3.8  km  and  10  km,  and  highly 
unstable  near  2  km.  This  instability  could  enhance  the 
gravity  wave  at  low  altitudes,  and  the  resultant  stability 
from  3.8  to  10  km  is  not  enough  to  damp  it. 

One  can  also  compare  the  two  soundings  generated  by  the 
MENTOR  to  see  if  the  differences  can  be  attributed  to 
meteorological  phenomena.  In  comparing  the  horizontal  speed 
along  with  its  zonal  and  meridional  components,  one  can  see 
evidence  of  a  horizontally-propagating  mountain  lee  wave. 
It  appears  that  the  zonal  component  was  perturbed  in  the 
first  profile  between  3-5  km  and  9-11  km.  These  altitudes 
are  where  the  zonal  component  becomes  easterly.  The 
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Fig.  37.  A  comparison  of  vertical  velocities  determined 
by  the  MENTOR  18:45:10-18:55:10  data  and  the  AFGWC 
HIRAS  model  for  Denver. 
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Fig.  38.  Plot  of  baroclinicity  vs 
MENTOR  18:45:10-18:55:10  data. 
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meridional  component  also  shows  evidence  of  this  wave  at  the 
same  altitudes  in  the  first  profile.  This  is  evident  by  the 
shift  to  a  more  northerly  component  at  these  altitudes.  The 
strong  directional  shear  also  points  to  a  wave.  It  appears 
that  the  wave  was  propagating  vertically  also;  as  the  second 
profile  of  vertical  velocity  shows.  Coincidentally,  the 
amplitude  peaks  are  at  3  km  and  11  km. 

A  quick  error  analysis  of  the  two  profiles  showed  that 
the  average  standard  deviations  of  VH,  u,  and  v  were  4.35 
m/s,  4.64  m/s,  and  4.14  m/s,  respectively,  while  that  of  a 
was  9.83*.  These  are  well  within  acceptable  meteorological 
standards,  especially  considering  the  wave  activity.  The 
averge  standard  deviation  of  w  was  0.14  m/s,  which  is 
exceptional.  This  further  verifies  that  MENTOR  provides 
extremely  accurate  vertical  velocity  measurements.  The 
average  standard  deviation  of  S  was  0.85*.  While  there  is  no 
standard  by  which  to  compare,  the  fact  that  the  vertical  and 
horizontal  components  compared  well  indicate  MENTOR  also 
provides  accurate  baroclinicity  measurments. 

Unfortunately,  the  lack  of  continuous  data  and  the 
amount  of  continuous  data  needed  to  generate  wind  profiles 
makes  it  impossible  to  show  plots  of  winds  with  time.  These 
would  provide  the  chance  to  see  changing  winds  due  to  various 
meteorological  phenomena  as  well  as  atmospheric  wave  motions. 
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CHAPTER  V 

CONCLUSIONS  AND  RECOMMENDATIONS 

One  can  conclude  from  the  results  that  the  MENTOR  does 
indeed  provide  reasonable  wind  profiles  of  the  lower 
atmosphere.  The  results  are  not  as  accurate  as  one  would 
like,  but  this  radar  is  new,  and  like  anything  new,  it  must 
be  tested  and  developed  to  further  refine  it  and  to  make  it 
useful.  It  was  not  originally  designed  to  measure  lower 
atmospheric  winds,  but  to  use  meteor  echoes  to  measure  upper 
atmospheric  winds.  All  the  profiles  seem  to  indicate  that 
the  MENTOR  is  not  subject  to  the  finite  range  volume  effect 
described  in  Chapter  3.  If  it  were,  there  would  have  been 
large  horizontal  wind  shears  measured  over  a  very  small 
altitude  range,  on  the  order  of  40  meters  per  second  per 
kilometer  (Fukao  et  al.,  1988).  However,  this  effect  was 
first  discovered  on  an  MST  radar.  The  Platteville  ST  radar 
did  not  appear  to  be  affected  by  it  either.  More  than 
likely,  the  conditions  were  not  right  for  this  effect  to 
occur. 

The  problem  of  the  distinct  power  drop-off  at  more  than 
2°  off-zenith  also  hinders  the  calculation  of  horizontal 
velocities.  A  small  decrease  in  frequency  would  provide  a 
wider  beam  and  could  provide  better  horizontal  velocity 


measurements . 
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This  decrease  in  frequency  could  allow  the  collection 
of  more  useful  points,  which  would  cut  down  the  number  of 
consecutive  soundings  needed  to  provide  a  profile  as  well  as 
increase  the  vertical  resolution.  The  frequency  decrease  may 
sacrifice  some  of  the  lower  altitude  information,  but  the 
increased  resolution  could  account  for  most  of  that. 

Another  area  of  potential  is  the  development  of  a 
desmearing  routine,  which  would  take  into  account  the 
movement  of  scatterers  across  different  Fourier  windows  as 
they  move  through  the  radar  beam.  This  would  improve  the 
accuracy  of  the  winds  that  are  calculated. 

There  is  also  the  problem  of  the  missing  y-component  on 
one  of  the  antennas.  Although  filtering  of  the  data  probably 
removed  most  of  the  bad  points,  fixing  this  problem  could 
also  improve  the  accuracy. 

I  believe  that  the  calculation  of  baroclinicity  could 
be  useful  to  determine  the  atmospheric  stability  at  various 
times,  and  could  be  a  useful  input  parameter  to  baroclinic 
numerical  weather  prediction  models. 

With  the  MENTOR  in  its  present  configuration,  I  found 
that  it's  best  to  use  a  10-minute  sounding,  a  linear 
polarization  filter  in  conjunction  with  doing  an  FFT  on  the 
time  domain  data,  filtering  points  using  a  least-squares 
polynomial  fit  to  a  plot  of  radial  zenith  angle  versus  radial 
velocity  calculated  by  coherently  averaging  over  one  minute, 
and  finally  filtering  out  all  points  outside  ±2°  off-zenith. 
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A  zenith  angle  window  of  24°  worked  best,  and  anything  less 
than  that  just  doesn't  allow  in  enough  points  to  plot  an 
adequate  profile.  The  inclusion  of  extremely  fast  data 
processing  hardware  and  software  will  make  it  possible  for 
MENTOR  to  generate  these  profiles  quickly,  thus  making  them 
much  more  current  and  therefore  more  useful  to  those  that 
require  real-time  wind  information.  The  current  rawinsonde 
network  only  provides  this  information  every  12  hours;  the 
data  is  old  30  minutes  after  it's  gathered. 

As  a  final  note,  if  all  these  suggestions  do  indeed 
produce  good  results,  then  the  MENTOR  can  then  be  included 
with  the  ST  and  SAD  methods  for  measuring  atmospheric  winds. 
There  are  also  possibilities  for  detecting  wave  motions  as 
well.  Best  of  all,  the  MENTOR  is  much  smaller,  requires  less 
power,  and  is  more  transportable  than  either  of  the  others. 
It  could  become  a  cost-effective  alternative  to  the  expensive 
rawinsonde  network  and  the  other  radars. 
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Fig.  4.  Complex  Fourier  spectra  for  the  two  antennas 
in  Figure  3,  showing  the  signature  of  a  single  moving 
target.  (From  Adams,  et  al.,  1985.) 
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